QUASI-PERIODIC SOLUTIONS FOR 
MATRIX NONLINEAR SCHRODINGER EQUATIONS^ 
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ly^ , The Adler-Kostant-Symes theorem yields isospectral hamiltonian flows on the dual 

0+* of a Lie subalgebra £|+ of a loop algebra g. A general approach relating the method of 
integration of Krichever, Novikov and Dubrovin to such flows is used to obtain finite-gap 
^..^ solutions of matrix Nonlinear Schrodinger Equations in terms of quotients of ^-functions. 
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In a recent paper [AHHl] the method of integrating nonlinear partial differential 
equations due to Krichever, Novikov and Dubrovin (see, for example, [KN, D]) was 



D ! adapted to computing finite-gap solutions for PDE arising as integrability conditions for 

'^ ' Lax equations arising in the framework of the Adler-Kostant-Symes (AKS) theorem. 

Consider a loop algebra g = 0"*" © 0~, split into the direct sum of the subalgebra 
^ ■ 0"*" of loops ^(A), which, viewed as maps X : S'^ \-^ g into a Lie algebra g extend 

holomorphically inside the unit circle in the complex A-plane and g", the subalgebra 
of loops extending holomorphically outside the unit circle, normalized by the condition 
X(oo) = 0. The Lax equations considered in [AHHl], such as (1.5a,b), determine 
completely integrable Hamiltonian systems. They follow from the AKS theorem and 
are are defined on the dual space g"""* of the subalgebra g"*". 

In an earlier paper [AHP], flows in a finite dimensional subspace consisting of 
rational elements in g+* were related to isospectral rank-r deformations of a fixed 
N X N matrix. In [HW] the coupled nonlinear Schrodinger equations corresponding 
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2 M.A. WISSE 

to the classical Lie algebras (cf. [FK]), which will henceforth be called "matrix NLS 
equations" , were related to the [AHP] framework. The general form of these equations 
is 

V^qt = qxx + qpq 

-\^Pt =Pxx+pqp, 

where q,p'^ G C"^^. 

Let sl{r, C), r = a + b he the Lie algebra of loops in 3l{r, C). Let sl{r, C)+ be 
the subalgebra of loops extending holomorphically inside the unit circle. Using the 
ad-invariant inner product 

^^y.^^I tjimnxn,, x,Yei(r,C) (1.2) 

we may henceforth identify the regular part of si{r,C)* with s[(r, C). In the same 
way sl(r, C)"*"* is identified with 5l(r, C)(^, the subalgebra of loops extending holo- 
morphically outside the unit circle. No notational distinction will be made between 
elements of sl(r,C)* (resp. sl(r, C)+*) and elements of 5l(r, C) (resp. 5l(r, C)(^). Let 
^_i_ = I{sl{r, C)*)~[/ ^-,^+* be the ring of Ad*-invariant functions on sl{r, C)* restricted 

to sl(r, C)"'"*. The AKS-theorem tells us that elements of JF_|_ commute in the Lie- 
Poison structure of sl(r, C)"'"*. Furthermore, the Hamiltonian equations for functions 
in JF_|_ are given by Lax pairs. The Hamiltonians for the matrix nonlinear Schrodinger 
equation are given by 

$.(X) = ltr((-^X(A))J (L3a) 



2 \\\' 



MX) = ^tr{{^XiX)) ) (L3b) 



for the X and t flows respectively, where X G sl{r, C)* and a(A) = nr=i(-^ " ^i)- "^^^ 
ai are constants chosen inside the unit circle. (N.B. This restriction is inessential; we 
can always choose our circle S^ to have a larger radius, so as to enclose arbitrarily 
chosen a^'s). Let 

where the A/^ G si{r, C) are of flxed rank ki < r — 1. Hamilton's equations for (1.3a,b) 
for elements of the form (1.4) are given by 

^Af{X) = [d^,{Af) + ,An (L5a) 

j^Af{X) = [dM^) + ,A^ (L5b) 
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where the "+" subscript means projection to sl(r, 0) + . The co-adjoint orbit Oj\j- 
through N preserves the pole structure of A/'(A) and the rank of Ni and hence these 
are preserved under any Hamiltonian flow. This orbit is therefore finite dimensional. 
The Lax equations for 



C{\) = ^U{X) = A"-i£o + >r~''Ci + ■■■ + U-x (1.6) 

A 



are given by 



^£(A) = [A£o + A,/:] (1.7a) 

^£(A) = [a2/:o + A/:i + £2, /:] (i.rb) 

These equations may also be viewed as due to the AKS theorem on the orbit On 
through £, with Hamiltonians 

$.(X) = ^tr((A2-X(A))o) (1.8a) 

$t(X) = ltr((A3-X(A))o) (1.8b) 

Equations (1.1) are the integrability conditions for (1.7a,b) if the leading terms of £(A) 
are 

. V^^ (hia A .^ „ . 

qp -q^ 



£2 = V-l n^ ^M (1.9c) 

\Px -pqj 

The underlying constraints on Cq, £1 and £2 are invariant under the fiows of all Hamil- 
tonians in JF^. 

Section 2 contains a brief summary of the main results of [AHHl]. The first part 
is devoted to the properties of the invariant spectral curve arising from Lax equations 
of the form 

£. = [P(£(A),A-i) + ,£(A)] (1.10) 

where C{X) is a matrix polynomial of the form (1.6), r is the fiow parameter and 
P{z,X~^) is a polynomial in z and A~^. The second part of section 2 summarizes 
the integration method for equations of type (1.10) yielding solutions in terms of 6- 
functions. Section 3.1 discusses the singularities of the spectral curve underlying equa- 
tions (1.7a,b). Solutions to (1.1) in terms of quotients of ^-functions are obtained in 
section 3.2. 
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Althogh only the generic case of a x b dimensional rectangular matrices q,p^ is 
treated here, we emphasize that finite-gap solutions for matrix NLS equations cor- 
responding to the various reductions by involutive automorphisms of the hermitian 
symmetric Lie algebra (5l(r, C),sl(a, C) © sl(6, C) © C) (see [FK, HW]) can be ob- 
tained by imposing appropriate initial conditions on the matrix polynomial (1.6). (See 
also [AHHl] for a general approach to reductions to sublagebras of 5l(r, C)). 

2. Preliminaries 

In this section we summarize the basic results of [AHHl]. 

2. 1 The spectral curve for Eows in a finite dimensional rational co-adjoint orbit 

This subsection is devoted to the specific properties of the spectral curve Sq given 
by 

V{z,X) = det{C{X)-zI)=0. (2.1) 

which is invariant under fiows of the Lax equation (1.10). 

Let 0{i) denote the i-th power of the hyperplane bundle 0{1) over P^(C). The 
technique of [AHHl] consists in compactifying Sq by embedding it into the surface 
T = 0{n — 1). On P^(C) consider standard coordinate charts (Vo,A) and (Vi,A) 
over Vo = pi(C)\{cx)}, Vi = pi(C)\{0}, with A = X'^ on Vo n Vi. The line bundle 
TV : 0{n — 1) ^ P^(C) is covered by the coordinate charts Ui = 7v~^{Vi), z = 1, 2 with 
coordinates (A, z) on Uq and (A, z) on t/i, where z = zX~^'^~^') on Uq fl Ui. Expanding 
(2.1) in z and A shows how Sq may be embedded into Uq- Changing to {X, z) over 
Uq n Ui extends the embedding of Sq into Ui and determines a compact curve S on 
the surface T that coincides with Sq over Uq- Using the adjunction formula [GH, p. 
471] the virtual genus 5' of 5 is easily computed to be 

g = ^{r-l){r{n-l)-2). (2.2) 

Since C is given by (1.6), the curve has singularities at the points (A, z) = (a^, 0), i — 
1, . . . ,n if ki =rank(A/i) < r — 1. Generically we have an ordinary (r — ki)-fold inter- 
section over the points A = a^. This is the case when the matrices A/^, z = 1, . . . , n are 
diagonalizable. Blowing up S at these points we obtain the (possibly singular) curve 
iS, with virtual genus 

^,,_f; C--^.)C'-^.-') , (2,3) 

2.2 Periodic Solutions in terms of 6 -functions for AKS-fiows 

In what follows, it is assumed that the partially desigularized curve S does not 
have other singularities. In the case of equation (1.1), however, there are generically two 
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additional singular points at A = oo. Section 3 shows how solutions can be computed 
considering these singular points. 

Let O-rii) = 7r*(9(z) denote the puUback of 0{i) to T. The flow of a matrix 
polynomial (1.6) gives rise to a flow of sheaves E over T, deflned by the exact sequence 

— ^OTi-n + l)®"^^^^"^ Of"^ ^E — .0. (2.4) 

where Ot is the sheaf of holomorphic functions on T, identifled with the trivial line 
bundle. By a standard abuse of notation we make no distinction between a line bundle 
and its sheaf of sections. 

If S is nonsingular and irreducible with genus g, E is a, line bundle of degree 
g + r — 1 over S. If S is singular and irreducible, E is a pushdown by some (possibly 
partial) desingularization map \1/ : iS — > 5 of a line bundle of degree g + r ~ 1 on S. 

Conversely, a linear flow of line bundles Ej- over S gives rise to a flow of matricial 
polynomials (1.6), governed by a Lax equation (1.10). 

Let C{X; 0) be an initial value polynomial deflning a spectral curve S. Let S be 
the desingularization of S with genus g. Let Eq be the initial value line bundle over S 
deflned by (2.4). The degree of Eq is generically g + r — 1. If Et- is a line bundle of 
degree g + r — 1 undergoing linear flow and having Eq as initial value, the tensor product 
Fr = Eq ^ Er is a line bundle of degree zero and hence has a transition function goi 
from Uq to Ui given by an exponential, 

goi{z, A) = exp(r;u(^, A"^)), (2.5) 

where fi is a, polynomial in z and A~^. Fix a basis {i/jI, . . . , '0^} of sections H'^iS, Er) 
(viewed as functions over a neighborhood of A = oo by flxing a local trivialization) 
normalized by the condition 

V';(oo,) = 5i, Vr, Vi,j = l,...,r. (2.6) 

where {ooi, . . . , oo^.} G S are the r points over A = oo. For pj in a neighborhood of ooj 
and A = 7r{pj) its projection to P^(C), deflne the matrix V't(A) by 

{Mm = i^UPj)- (2.7) 

Let z{pj), j = 1, . . . , r be the (generically distinct) eigenvalues of £(A; 0). The matrix 
polynomial C{X) given by 



£(A;r)=^-i 



/4pi) \ 



\ Z{pr) j 



tpr (2.8) 
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then satisfies the Lax equation 

^C{X; t) = [^(£(A; r), A-i) + , £(A; r)] (2.9) 

where the "+" subscript means projection to positive powers of A. 

Remark. In the AKS framework, fi{z,X~^) is given by d(f)r{z), where (pr is the 
Hamiltonian defining the r-fiow. 

To obtain a representation of the sections t/j^, j = 1, . . . ,r in terms of 6 functions, 
we first describe exphcitly the sections of Eq. Let A be an effective divisor representing 
Eg. For instance, picking a fixed vector v G C^ we can take A to be the sum of points 
A^ = {Xj,Zj)j=i^,„^g+r-i in5 such that {C{\j)-ZjI))adjV = 0, but {C{Xj)-ZjI)adj 7^ 0. 
Here the subscript denotes the classical adjoint matrix and the vector v G C^ represents 
an element of H^{S, Eq). Let D^o = C)Oi + ■ ■ ■ + oo^ be the divisor consisting of the 
r points over A = oo. Due to the Riemann-Roch theorem, assuming A — Doo + ooi 
is in general position there is exactly one meromorphic function i/jq on S such that 
(^o) ^ -^ + -^oo - ooi and 

V^o(oo,)=<5]. (2.10) 

Such a function is also a section of Eq. The r sections {ipo, ■ ■ ■ , V'o) form a basis of 
H^{S,Eo). Following [AHHl] such sections may be expressed in terms of quotients 
of 6'-functions defined on the Jacobi variety J' of S. 

Choose a basis of cycles (ai, . . . ,ag,bi, . . . ,bg) onS with the intersection property 
tti ■ ttj = 0, bi ■ bj = and ai ■ bj = 6ij, i,j = 1, . . . , g. Let {a;i , . . . , ujg} be a basis of 
holomorphic differentials on S, normalized by the condition that §^.uJi = <^o-, hj — 
1, . . . ,g. Define the Abel map A : S x S -^ J' by 

A{x,y)= ( [ ui,..., [ Ug). (2.11) 



The corresponding 6 function is defined by the matrix of 6-periods of {uji, . . . ,Ug}. 
Define Q} + ■ ■ ■ + Qf = {%) - D^ + ooi + A, with A = A^ + ■ ■ ■ + A^+'^-i. There 
exists e G C^ such that 

9{e) = 

^(e + A(oo„y))^0, Vj = l,...,r 
e{e + A{A^,y))^0, Vj = 1, . . . ,^ + r - 1. 
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Fix a base point p G 5 and define constants 

i=i i=i 

r 

(3 = J2Mp,oOj) (3, = A{p, oo,) (2.13) 

^. = J2Mp,Qi)- 

Let 5 be the Riemann constant. With these definitions the functions 

e{A{p,y) + 5-a'-a" + P-P,)Uj^^0{A{oo,,y) + e) 

F (y) = i — (2.14) 

9{A{p,y) + 6-a')U'jZlO{A{A-s+^,y) + e) 

define sections of Eq with the correct zeros and poles. The normahsation condition 
(2.10) is obtained by setting 

FHy) 
%iy) = ^-^. (2.15) 

To compute sections of the time dependant hne bundle E^- = Eq ® Fj- we recall 
that Ft is given by the exponential transition function (2.5). Thus sections of E^ 
can be represented as functions with zeros at D^o — oo^, poles at A and exponential 
singularities at -Dq = Oi + ■ ■ ■ + 0^-, the r points over A = 0. Such functions are called 
r-point Baker-Akhiezer functions. 

For an explicit representaion of ip^. let /I be a differential of the second kind on S 
with zero a-cycle integrals and the same polar part as d^i. Let U G C^ be given by 

Define 

K(y) = exp [r j^ ,.j ^(A(p.,) + 6-,.) '"^' 

then the functions _ 

/i^(oo,)F'(ooO 
define a basis of H^{S,Er) satisfying properties (2.6). 
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3. The Matrix Nonlinear Schrodinger Equation 

The techniques of section 2.2 also apply to the matrix NLS equation once the 
corresponding spectral curve is desingularized at A = oo. 

3.1 The Spectral Curve for the Matrix NLS Equation 

For the matrix NLS equation the leading terms of C{X) are given by (1.9a-c). 
In order to study the behaviour near A = oo of the corresponding spectral curve S, 
given by (2.1) embedded into T, we switch to coordinates (A, 5) which leads to the 
representation 

P(A, z) = det(£(A) - zl) = 0. (3.1) 

with £(A) given by 

£(A) =Co + A£i + ~X^C2 + ■■■ + A"-^/:„-i. (3.2) 

Expanding (3.1) in coordinates {X, z' = —{z + V~l^^)) around {X, z') = (0,0) gives, 
under invariant genericity conditions on £3, 

6 a 

l[{z' + Prf, + 0(A^)) Y[{z' + e, + 0(A2)) = (3.3) 

where dj,ei are constants. This shows that S has a fe-fold intersection of order 3 at 
(A = 0, 5 = — V— 1^^)- Blowing up three times at this point reduces the virtual genus 
of S by — 2 Similar things happen at the (generically) a-fold point of order 3 in 
(A = 0, 5 = V— 1^^)- Blowing up three times at this point reduces the (virtual) genus 

1 3a(a— 1) 

by \ > . 



Call S the curve obtained after desingularization at the n + 2 points (a^, 0) 



i=l, 



(00, — a/— T^r), (00, a/— T— ^), in coordinates (A, z). Invariant genericity conditions on 



L may be chosen such that S does not have other singularities. The arithmetic genus 
^ of 5 is therefore 

~g=]^Ur- l)(r(n - 1) - 2) - f ^(r - k,)(j- - A;, - 1) j - 36(6 - 1) - 3a(a - 1) 

(3.4) 
3.2 Solutions for the Matrix NLS Equation in terms of 9 -Functions 



In order to use formulas (2.16-18) we need an explicit representation of the poly- 
nomials ii{z, X~^) and v{z, X~^) defining transition functions qqi from Uq to Ui for the 
degree zero line bundle F^ t by 

5^01 (A, z) = exp{xn + tv) . (3.5) 
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As mentioned before, in the AKS framework n and u are given by the differentials of 
the AKS hamiltonians (1.8a,b). An easy computation gives 






A 



n-2 



Xn-3 



(3.6a) 
(3.6b) 



The (x, t)-dependance for the matrix polynomial (3.2) is computed as in (2.8). 
The components of £i carry the solution for the matrix nonlinear Schrodinger equation 
(1.1). It is recovered by differentiating (3.2) with respect to A at A = 0. The solutions 
q, p to (1.1) are given by 



dX 

p){x,t)= -^ 
d\ 



C{\) 



j+a 



A=0 



£(A) 



A=0 



j+a 



(3.7a) 
(3.7b) 



i = 1, . . . , a, J = 1, . . . , fe. A short computation yields 



qj {x, t) 
p){x,t) 



-l(^^'+")'(cx) 
^(V^^)'(oo 



j+o.) 



(3.8a) 
(3.8b) 



where the prime designates derivation with respect to A. Let /x, v be normalized dif- 
ferentials of the second kind on S with same polar part as li/x, dv respectively. Set 






(3.9a) 
(3.9b) 



Fix a point p E S and substitute formulas (2.14,2.17) into (2.18). Use (3.8a,b) to get 



qi{x,t) = Ki+''exp{eix + f^t) 



e{A{p, ooi) + xU + tV + S- 7j+c 



eiA{p, OOj+a) +xU + tV + S- -fj+a] 



p]{x,t) = -i^j+,exp(-efa: - //t) 



e{A{p, ooj+g) +xU + tV + S--fi 
e{A{p, oo^) + xU + tV + 5- 70 



(3.10a) 
(3.10b) 
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where the constants Kj', el and ff are given by 

-e{A{p, oofc) + 5 - ik)e{A{p, ooz) + 5 - a' - a" + /? - Hi) 



Kt = V-1- 



F''{ook)e{A{p,ooi) + 5-a') (g^^^) 

e{A{ooi,y) + e) T.m^k,i 0{A{oom, ooi) + e) 



_d_ 

dy 



y=ooi 



OOi 



n:rii^(^(A^+"^,oo,)+e) 

ei= I fl (3.11b) 

JoOj + a 

n= ^ (3.11c) 

where the integrals are evaluated in the same polygonization as that defining the Abel 
map A. 

The functions (3.10a,b) are the desired periodic solutions to the matrix nonlinear 
Schrodinger equation (1.1). 

Comments 

The solutions for the matrix NLS equations obtained above could also have been 
obtained by applying the Hamiltonian integration method based upon computation 
of the Liouville- Arnold generating function (see [AHH2,3]) determining a canonical 
transformation to linearizing coordinates. That method explicitly shows how nonlin- 
ear Lax equations (1.7a,b) linearize on the Jacobi variety associated to the invariant 
spectral curve and makes the evaluation of many of the constants more explicit. How- 
ever, use of the Krichever-Novikov-Dubrovin method avoids problems associated with 
"missing Darboux coordinates" due to the singular spectral data at A = oo implied by 
the sructure of the leading terms in C given by (1.9a-c). In the Hamiltonian approach 
this requires the imposition of other symplectic constraints defining a reduced family 
of spectral data involving curves of the same genus as the desingularized ones. The 
main shortcoming of the present approach lies in the difficulty on explicitly computing 
the constants 7^ and K\ in view of the complicated implicit way they arise through 
formulas (2.13), (3.9a,b) and (3.11). 
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